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ABSTRACT 

The multiwavelength observation of the nearby radio galaxy M 87 provides a unique oppor- 
tunity to study in detail processes occurring in Active Galactic Nuclei from radio waves to TeV 
7-rays. Here we report the detection of 7-ray emission above 250 GeV from M 87 in spring 
2007 with the VERITAS atmospheric Cherenkov telescope array and discuss its correlation with 
the X-ray emission. The 7-ray emission is measured to be point-like with an intrinsic source 
radius less than 4.5 arcmin. The differential energy spectrum is fitted well by a power-law 
function: d$/dE=(7.4± 1.3 stctt ± 1.5 sys )x (E/TeV)^ 2 - 31±0 - 17st " t±0 - 2 ^)fO- 9 m- 2 s- 1 TcV- 1 . We 
show strong evidence for a year-scale correlation between the 7-ray flux reported by TeV experi- 
ments and the X-ray emission measured by the ASM/RXTE observatory, and discuss the possible 
short-time-scale variability. These results imply that the 7-ray emission from M 87 is more likely 
associated with the core of the galaxy than with other bright X-ray features in the jet. 

Subject headings: Galaxy: individual(M 87, Virgo A, NGC 4486) - gamma rays: observation - VERITAS 
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Introduction 

The giant elliptical galaxy M87 is a nearby 
-16 Mpc) powerful FR I (jFanaroff fc Rilev 



19741 ) radio galaxy (Virgo A) which lies near 
the center of the Virgo cluster. Its core is 
an Active Galactic Nucleus (AGN) powered 
by a sup ermassive black hole o f (3.2 ± 0.9) x 
10 9 M Q (jMacchetto et all Il997l). em itting the 
first-observed plasma jet ( Curtis! 19181 ). M 87 has 
been observed over a broad range of energies from 
radio waves to TeV 7-rays. Its jet is resolved in 
radio, optical and X-ray regimes and shows sim- 
ilar morphologies at all these wavelengths. M87 
is the first non-blazar AGN observed to emit TeV 
7-rays and it provides valuable insight into the ac- 
celeration of high-energy particles in astrophysical 
jets. 

The first detection of M 87 in the TeV regime 



was r eported by the HEGRA collaboration (jAharonian et 
2003) with a statistical significance of 4.1 stan- 
dard deviations derived from observations made 
during 1998-99. This de tection was confirmed b y 
the HESS collaboration (jAharonian et alJlioOGAf ). 
which also reported a year-scale flux variability 
and strong indications of rapid variability (2-day 
scale) during a high state of 7-ray activity in 2005. 
The imaging atmospheric Cherenkov technique 
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provides insufficient angular resolution (few ar- 
cminutes) to resolve the M 87 emission region, but 
the day-scale variability claimed by HESS sug- 
gests a very small region, most likely close to the 
core. The confirmation of such short-time-scale 
variability is of paramount importance because of 
its implications regarding the source dimensions. 

The rapid variability of M 87 7-ray emission 
would also constitute an important connection 
with the most common extragalactic sources de- 
tected in the TeV regime. These are all gener- 
ally considered to be BL Lacertae (BL Lac) ob- 
jects with the jet pointing towards the observer 
(blazar). The TeV emission of M87 is now gener- 
ally inter preted in the frame of th e unified scheme 
of AGN (|Urrv fc Padovanilll995lh According to 
this scheme, FR I radio galaxies are of the same 
nature as BL Lac objects, but with their jet not 
ointing along the line of sight. Typical mod- 
is assume a M87 jet misalign ment around 30°- 
35° (jBicknell fc Begelmanlll996l ). However, super- 
luminal motion has been observed i n the compact 



first knot, HST-1, in th e optical (IBiretta et al.. 
Il999h and radio domains (jCheung et al.ll2007l ). To 



explain this requires a jet orientation closer to the 
line of sight (within 19°) at least at this location in 
the jet. The HST-1 emission also shows "blazar- 
like" behavior with very strong month-scale vari- 
ability in the radio , optical and X-ray wavelengths 
(jHarris et alJ 120071 ). much stronger than the vari- 
ability observed in the nucleus or in the other 
knots. 

The non-thermal emission of the jet is un- 
derstood as synchrotron radiation from high- 
energy electrons, and the TeV 7-rays could re- 
sult from inverse-Compton scattering by the same 
electron population. A time correlation be- 
tween the X-ray synchrotron radiation and the 
TeV emission is then expected. Models explain 
the TeV emission as originating from electrons 
either in the inner jet close to the core and 
the jet emission base ( TeV b la zar-type model) 
Georganopoulos et al.l (120051) . iGhisellini et al 



2005). lLenain et alJ (|2007h . lTavecchio fc Ghisellin 
(12008 )), or in the large-s c ale je t (jStawarz et al.l 
(|2003l ). iHonda fc Honda I (120071 )). On interme- 
diate scale, it was recently proposed the TeV 
emission comes from the peculiar knot HST-1 
( Stawarz et al. 20061 ). Since the year-scale vari- 



ability is established, the core region and the knot 
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HST-1 are more likely the source of the TeV emis- 
sion than the large-scale jet. 

The origin of the TeV emission could also not be 
in the jet but in the vicinity of the central super- 
massive black hole of M 87. The rotating magneto- 
sphere of the black hole could accelerate particles 
(pulsar-type model) producing electromagnetic 
cascades a nd then TeV 7-rays by inverse - Compton 
scatte ring ([Neronov fc Aharonian (|2007h .lR ieger fc 
This model could explain a day-scale 



tions from 55° to 71°. In spring 2007, VERITAS 
was still in its construction phase. The M87 ob- 
servations were performed with an array of either 
three telescopes (94% of the data) or four tele- 
scopes (6%). The performance of th e array dur 



ing th is period is discussed elsewhere (jMaier et al 
In order to minimize systematic errors on 



the estimation of the flux and its variability, we do 



variability of the TeV emission. 

An alternati ve model involving p rotons has also 
been proposed (jReimer et al ]|2004h . In this model, 
the TeV 7-ray emission is dominated by neutral 
pion production, and by proton and muon syn- 
chrotron radiation in the highly magnetized envi- 
ronment of the jet formation region (Synchrotron- 
Proton-Blazar model). In this context, M87 
would be an efficient cosmic rays accelerator. It 
was even suggested that M 87 is the source of most 
ultra-high-energy cosm ic rays detected on Earth 



AharorBSt] use the data from the fourth telescope, which 
was not fully calibrated at this time. We processed 
all the data in the same way, using only three tele- 
scopes, even when data from the fourth telescope 
were available. 

Only good-quality data were considered for 
the analysis. Our data selection is based on 
weather conditions and raw trigger rate stabil- 
ity. About 90% of the data (44 h) pass our quality 
selection cuts. The selected data have been pro- 
cessed with several independent analysis packages 



(jBiermann et al.l l2001). However the SPB model 



predicts a quite soft spectrum in the TeV range, 
softer than what was recently measured by HESS 
and VERITAS. 

Another model, involving dark matter annihila- 
tion, has also predicted TeV 7-ray emission from 



M87 (|Baltz et al.l l2000h . But the variability of 



the flux almost excludes this scenario as the main 
contributer to the TeV emission. Dark matter is 
expected to be distributed over much larger dis- 
tance scales than those implied by the recently 
established year-scale variability. 

Studies of the variability time scale and of coin- 
cident multiwavclcngth observations are very im- 
portant for discriminating between these models. 

2. Observations and analysis 

The Very Energetic Radiation Imaging Tele- 
scope Array System (VERITAS) is a 7-ray ob- 
servatory located in southern Arizona. It is 
an array of four 12m-diameter imaging atmo- 
spheric Cherenkov telescopes, each with a 499- 
pixel photomul tiplier-tube camera covering a 3.5° 



field of view ([Holder et al.l l2006h . The 



energy 

range covered extends from 100 GeV (zenith ob- 
servation) up to tens of TeV. 

The observations of M 87 were carried out over 
51 hours between February to April 2007 at eleva- 



(jDaniel et al.l 120071 ). The results presented here 



come from one of these packages and the others 
yield consistent results. The analysis chain uses 
two-threshold image cleaning, stereo reconstruc- 
tion and cosmic-ray re jection based on standard 
image parameter cuts (jHillas 1 11985). The energy 
of each shower is reconstructed using a method 
based on the quantity of light measured by the 
telescopes and on the position of the shower-core 
ground impact. The differential energy spectrum 
is unfolded using the correction-factor method 
described in Chapter 11 of "Statistical Data 
Analysis" (jCowanll 19981 ) . 

All observations were taken in "wobble" mode, 
tracking M87 with a 0.5° offset north, south, east 
or west from the camera center. In this mode the 
cosmic-ray background rate is estimated using sev- 
eral "mirror" regions at the same distance from the 
camera center (0.5°) as the nominal source posi- 
tion. We compare the number of events with the 
reconstructed shower direction in these off-source 
regions (seven in this analysis) to the number of 
events in a region centred on the nominal source 
position. The radius of each region is 0.13°. Sta- 
tistical significances are calcul ated using the Li 
and Ma formula (equation 17 in lLi k MaT(|l983h ). 



3. VERITAS results 

We detected an excess of 259 ± 44 events from 
the direction of M87, corresponding to a statis- 
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tical significance of 5.9 standard deviations. The 
analysis energy threshold for the average elevation 
of the M87 observations (66°) is approximately 
250 GeV, estimated from Monte Carlo simulations. 
The time-averaged excess rate (0.10±0.027/min) 
corresponds to 1.9% of the Crab Nebula rate ob- 
served at similar elevations. 

Figure 1 shows the map of the excess signif- 
icance in the sky region of M87. The detected 
source (VER J1230+123) is consistent with a 
point-like source located inside the M87 galaxy. 
The position of the maximum of a 2D Gaussian 
fit of the excess is: R.A. 12 h 30 m 46 s ± 4 s stat ± 6 s sys , 
Decl. +12°23'21" ± 50' s ' tcri ± l'30^ s , compat- 
ible with the position of the M87 core (R.A. 
12 /l 30 m 49.4 s , Decl. +12°23'28"). Figure 2 shows 
the distribution of the square of the angle 9 be- 
tween M 87 and the reconstructed shower direc- 
tion. The shape of the excess is compatible with 
a point-like source (dashed line). The upper limit 
(99%) of the intrinsic source extension has a 4.5' 
radius. 

The differential energy spectrum has been mea- 
sured from 250 GeV to 8TeV. Table 1 gives the 
differential flux and statistical error at different 
energies. In order to keep the reconstructed en- 
ergy bias below 10%, the flux of the lowest en- 
ergy bin (250 GeV) has been estimated using only 
data taken at elevation above 62° (34.8 h of ob- 
servation). The M87 spectrum obtained is consis- 
tent with a power-law spectrum: d$/dE=<I>o • 
(_E/TeV) r with a flux normalization constant 



4>n 



(7.4 ± 1.3 stai ± 1.5 fl „ fl )10- 



TeV~ 



and a spectral index T — — 2.31 ± 0.17 stat ± 0.2 sys . 
The x 2 /d.o.f. of the fit is 3.0/4. Figure 3 shows 
this spectrum and its power-law fit and compares 
them with the pre vious spectra measured by HESS 



in 2004 and 2005 (jAharonian et alJ 12006 At ). Data 



point s from a recen t re-analysis of the HEGRA 
data ( Gotting 20061 ) are also shown. The results 



of this HEGRA re-analysis will be used hereafter. 
Correlation between the spectral index and the 



as Mrk501 (Diannati et a 


.1 1999) and 


(Krennrich et all (20021). 


Albert et all 



(2007)). 



Unfortunately, the 7-ray flux from M87 being so 
close to the sensitivity of the present generation of 
atmospheric Cherenkov telescopes, the uncertain- 
ties on the spectral index do not permit to make 
any statement on a possible similar behavior. 
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Fig. 1. — Significance map of the region surround- 
ing M87 measured with VERITAS (in units of 
standard deviation). Circles show the positions 
of the brightest galaxies of the Virgo cluster. 
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Fig. 2. — The data points with error bars show the 
2 distribution from M 87 and the solid histogram 
that from off-source directions. The dashed line is 
the 2 distribution expected for a point-like source 
with a background level represented by the dotted 
line. 
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Fig. 3. — Spectral energy distribution of M87 
measured with VERITAS in 2007 and with other 
experiments during the previous years. The 
dashed line shows the power-law fit of the VERI- 
TAS data (see text). 



The 7-ray flux shows significant variability from 
one year to another. The first panel of Fig- 
ure 5 shows the integral 7-ray flux reported by 
the HEGRA, HESS and VERITAS collaborations. 
For case of comparison, the VERITAS flux was 
scaled to the energy threshold of 730 GeV com- 
monly used by HEGRA and HESS. The integral 
flux above 730 GcV is estimated according to our 
spectral analysis result (power-law spectrum with 
a spectral index r=-2.31): 
$>730GcV=(8.5±1.5 stat ±1.7 sys )10- 9 m- 2 s- 1 . 

The TeV 7-ray lightcurve consists of data 
points from different experiments. The pre- 
liminary result for the Crab Nebula, standard 
candle of the TeV astrono my, observ ed with 



VERITAS in spring 2007 (jCelik et all 120071 ) 
is in good agreement with the fl ux reported 
by HEGRA (lAharonian et all 12004) and HESS 



(jAharonian et alj I2006BI ) . Thus, the differences 



observed with different detectors comes from real 
variations of the M 87 flux. 

Figure 4 shows the M87 7-ray flux night-by- 
night during the three months of observation with 
VERITAS. No significant short-time-scale vari- 
ability is observed. The constant flux fit % 2 per 
degree of freedom (d.o.f.) is 24.3/22. The maxi- 
mum deviation from the average (MJD 54173) is 
only 2.2 standard deviations. 



ri5o = 



>100 

0) 

(3 



90 100 110 
Time(MJD-54100) 



Fig. 4. — Nightly average (circles) and global av- 
erage (dashed line) of the M87 7-ray flux above 
250 GeV measured with VERITAS as a function 
of the day of year 2007. 



4. X-ray / 7-ray correlation 

The All-Sky Monitor (ASM) on the Rossi X- 
ray Timing Explorer (RXTE) has been monitoring 
the M 87 e mission between 2 k eV and 10 keV since 
early 1996 ( Levine et al.lll996h . Measurements are 
performed with sequences of 90-second "dwells". 
The light-curve data used here are the quick-look 



Tabic 1: Differential flux of M87 measured with 
VERITAS in spring 2007 



Energy 


Differential flux 


Significance 


(TeV) 


(m- 2 s- 1 TeV- 1 ) 


(std dev.) 


0.25 


(1.36±0.60)10" Y 


2.26 


0.50 


(4.55±1.15)10~ 8 


3.97 


1.00 


(7.39±2.65)10- 9 


2.80 


2.00 


(1.92±0.68)10- 9 


2.82 


3.98 


(1.20±1.66) 10" 10 


0.73 


7.94 


(1.01±0.85)10- 10 


1.19 
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results provided by the ASM/RXTE team on their 
web sit43, dwell- by-dwell and as a daily average. 

In order to check the long-term variability of 
the ASM/RXTE signal, we sum the data by 
bins of 6 months using the standard practice of 
weighting measurements by 1/cx 2 , where a is the 
measurement uncertainty. The second panel of 
Figure 5 shows the average ASM/RXTE rate 
for the first six months of each year (i.e. the 
period coinciding with the observation of M87 
with the ground-based detectors) since 1996. The 
ASM/RXTE lightcurve has similar variations as 
the TeV lightcurve shown in the first panel. Fig- 
ure 6 shows the 7-ray flux, "3?>730GeV, plotted 
versus the ASM/RXTE rate, X rate . Taking the 
error bars, which correspond to the statistical 
error, into account, the linear correlation coeffi- 
cient is 0.78 ±0.11 (positive correlation confidence 
level > 99.99%). The most recent analysis of the 
HEGRA data has been used to calculate this cor- 
relation coefficient. 

A marginal correlation between the TeV 7- 
ray and the ASM/RXTE rate was previously 
suggested by the Whip ple lOm-telesc o pe col - 



laboration in 2000-01 (|LeBohec et all I2003AI ) 
and was discussed in the doctoral theses of 



M. iBeilicke I (120061) and of N. iGotting I (l2006h . 
Such a correlation between the TeV and hard 
X-ray emissions has bee n also detected in th e 
BL lac objects M r k501 (|Diannati et all (Il999l). 



Krawczvnski et all (|2002l )) and Mrk421 (jFossati et al 
2007). The correlation for these two objects is 
generally close to quadr atic but depends on th e 
energy range considered ( Katarzvnski et aLll2005f ). 

We estimate the value of both the variable and 
the steady components of the X-ray emission of 
M 87 by fitting the TeV correlation plot with a 
linear relation: 

$>730GcV=(88±23)-(V rate - (1.196±0.022))-10- 9 
and with a quadratic relation: 

*>730GcV=(338±180)-(V rate -(1.142±0.042)) 2 -10- 9 . 
In these two equations $>730GeV is expressed in 
m _2 s _1 and X rate in count/s. Both functions fit 
the data well with a x 2 /d.o.f.<l. They show that 
the X-ray emission consists of both a dominant 
steady emission (l.l-1.2count/s) and a variable 
fraction (<20%) strongly correlated with the 7- 
ray emission. 



The constant part of the X-ray emission is dom- 
inated by the thermal emission from the ~2keV 
gaseous atmosphere of the M87 galactic halo 



(jForman et al.l 120071 ). The variable part could 



result from synchrotron emission by high-energy 
particles with cooling time less than one year. The 
emission site for the variable component must have 
a compact size on the order of one lightyear or less 
(~5 mas at the M87 distance), or be moving at 
relativistic speed in a direction close to the line of 
sight, to explain the yearly variability recorded. 

The only known M 87 regions which have such 
characteristics are the core and the brightest knots 
of the jet. The Chandra X-ray observatory has 
been used to monitor the M 87 jet e mission in the 



0.2-6 keV energy range since 2000 (jHarris et al 



20061) with enough angular resolution to measure 



1 http://xte.mit.edu/asmlc/ASM.html 



separately the flux from the different features. 
Figure 5 shows the light curves of the two dom- 
inant X-ray features which are the core and HST- 
1. Since 2003, HST-1 has been flaring and it 
largely dominates the Chandra X-ray flux. How- 
ever the variation of the flux from this knot does 
not seem to correlate well with the variation of the 
ASM/RXTE signal. Between 2003 and 2005 the 
HST-1 X-ray emission increased by a factor of 20 
whereas the variable fraction of the ASM/RXTE 
rate changed by less than a factor of two. The 
same goes for the new flux measurement with 
VERITAS which clearly do not follow the X-ray 
lightcurve of HST-1. 

The next-brightest features of the jet observed 
with Chandra can explain only with difficulty the 
ASM/RXTE rate variation as their fluxes are sig- 
nificantly weaker than the HST-1 X-ray variation. 
Actually, the energy range of ASM/RXTE (2- 
lOkeV) is higher than that of Chandra (0.2-6 keV) 
and the dominant features of the ASM/RXTE en- 
ergy band may not be HST-1. According to the 
spectral analysis of Chan dra data taken in 2000 
(jPerlman k Wilson I [20051) . ;he core was the only 
bright X-ray feature with a harder spectrum than 
HST-1. Thus the core seems the best candidate to 
dominate HST-1 in hard X-rays and provide sub- 
stantial ASM/RXTE rate variation. Moreover the 
confidence level for a positive linear correlation be- 
tween the 6- month average rates of ASM/RXTE 
and Chandra is higher for the core (90%) than for 
HST-1 (65%). Since the variable component of 
the ASM/RXTE flux is likely dominated by X-ray 
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emission from the core, and there is a strong cor- 
relation between the TeV and ASM/RXTE fluxes, 
we conclude that the core is the most likely can- 
didate for the source of the TeV 7-ray emission. 

It is however surprising that the large HST-1 
flare in 2005 does not show in AMS/RXTE as this 
implies a very soft HST-1 spectrum. There is no 
evidence for such a soft spectrum in the Chan- 
dra data. Unfortunately, during the flaring pe- 
riod of HST-1, spectral analysis of the core and 
HST-1 seems impossible because the flux was so 
strong that the bulk of the events in the image 
are piled, and this destroys the spectral informa- 
tion. Furthermore the core is separated by only 
0.86" from HST-1, and its Chandra signal can be 
contaminated by t he HST-1 signal. As shown in 
Harris et al.l (J2007J), the small flare of the core in 
2005 can be almost smoothed out by subtracting 
5% of the HST-1 signal (dotted line in the third 
panel of the figure 5). This is probably indicative 
of the uncertainties on the core flux measurement 
with Chandra during this period. 

Motivated by the long-time-scale correlation 
between X rate and < i ) >730Gey, we searched for a 
similar correlation on a 5-day scale. Such short- 
time-scale correlation had previousl y been sug- 
gested from earlier Whipple 10m data ijLeBohec et al 
2003BI ). This study was also designed to im- 
prove our sensitivity to short-time-scale variations 
such as those reported by the HESS collabora- 
tion. Again, we express the VERITAS results as 
the integral 7-ray flux above 730 GeV. Figure 7 
shows $>730GeV versus X rate summed on con- 
tiguous and independent 5-day-scale bins. The 
confidence level for a positive correlation is only 
82%. The statistics are too poor for us to draw 
any conclusion as to a possible extension of the 
year-scale correlation (dash-dotted line) to the 5- 
day scale. Simultaneous observations of M 87 with 
VERITAS and ASM/RXTE will continue in the 
coming years, providing the opportunity to inves- 
tigate further such short-time-scale correlation. 

5. Conclusion 
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Fig. 5. — The first panel shows the integral 7- 
ray flux above 730 GeV reported by the HEGRA, 
HESS and VERITAS collaborations. Concerning 
the HEGRA flux, both the published result (light 
color) and the most recent re-analysis result (dark 
color) are shown (see text for references). The sec- 
ond panel shows the January-to-June average of 
the global X-ray emission of M 87 measured with 
ASM/RXTE in the range 2-10 keV. The last two 
panels show the M87 core and the knot HST-1 X- 
ray emission measured with Chandra in the range 
0.2-6 keV. The dotted line on the third panel shows 
the corrected core signal fro m which 5% o f the 
HST-1 signal were subtracted (jHarris Il2007h . 



The VERITAS collaboration observed M 87 in 
spring 2007 with an array of three telescopes. 
The analysis of the 44 hours of data provides a 
5.9 standard deviation detection of 7-ray emission 
above 250 GeV. The energy spectrum is fitted well 
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Fig. 6. — Correlation between M87 emission in 
X-ray (ASM/RXTE data: 2-10 keV) and in 7- 
ray (HEGRA, HESS, VERITAS data: >730GeV). 
The dash-dotted line is the linear correlation fit 
and the dotted line shows the quadratic correla- 
tion fit. For HEGRA, both the published flux 
(dashed lines) and the reanalysis result (solid lines) 
are shown. 
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Fig. 7 — M 87 X-ray emission (ASM/RXTE data: 
2-10 keV) plotted relative to the very-high-energy 
7-ray emission (VERITAS data: >730GeV) for 
contemporenous observations. Data are binned in 
a 5-day intervals. The dash-dotted line shows the 
linear correlation fit obtained for the year-scale 
correlation in figure 6. The grey contours are the 
sum of 2D Gaussian functions associated with each 
data point according to its uncertainty. 
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by a power-law function: d$/dE =(7.4 ± 1.3 stat ± 
1.5 sys )x (E/TeV)(- 2 - 31±0 - 17 — ±0-2. S a)i -9 m -2 s -i 
The comparison of our result with the spectra re- 
ported by the HESS and the HEGRA collabora- 
tions does not indicate a variability of the spectral 
index. 

In most TeV 7-ray emission models, a corre- 
lation between the 7-ray and X-ray flux is ex- 
pected and the relationship can be used to dis- 
tinguish between models. We have found strong 
evidence for such a correlation from M87. We 
have shown a year-scale correlation between the 
7-ray flux above 730 GeV, recorded by imaging 
atmospheric Cherenkov telescopes during the last 
ten years, and the X-ray flux in the energy range 
2-10 keV recorded with ASM/RXTE. Both linear 
and quadratic functions fit this correlation well. 
The poor correlation between the ASM/RXTE 
rate variations and the dominant X-ray features 
(HST-1) observed by Chandra in the 0.2-6 keV en- 
ergy band is surprising but, in the light of the cor- 
relation with the TeV flux, it suggests HST-1 is 
unlikely to be the main source of the 7-ray emis- 
sion. If the core indeed has a harder spectrum 
than HST-1, it could dominate in the ASM/RXTE 
energy range, and be the most likely site of the de- 
tected TeV emission. HESS observations of M 87 
in 2005 were suggestive of a 2-day-scale variability. 
The VERITAS observations in 2007 showed it to 
be in a lower state and, even with a short-time- 
scale X-ray correlation study, the statistics of this 
observation are insufficient to confirm such rapid 
variability. Further observations and multiwave- 
length campaigns on M87 at a range of activity 
states will prove important in tracking down the 
site of the high-energy emission, determining the 
particle acceleration mechanism and the relation 
to the blazar class of high-energy objects. 
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